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This paper develops a method of manipulating the 
squeezed atom state to generate a few-photon state whose 
phase or photon-number fluctuations are prescribed at our 
disposal. The squeezed atom state is a collective atomic state 
whose quantum fluctuations in population difference or col- 
lective dipole are smaller than those of the coherent atom 
state. It is shown that the squeezed atom state can be gener- 
ated by the interaction of atoms with a coherent state of the 
electromagnetic field, and that it can be used as a tunable 
source of squeezed radiation. A variety of squeezed states, 
including the photon-number squeezed state and the phase 
squeezed state, can be produced by manipulating the atomic 
state. This is owing to the fact that quantum-statistical in- 
formation of the atomic state is faithfully transferred to that 
of the photon state. Possible experimental situations to im- 
plement our theory are discussed. 
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I. INTRODUCTION 

Nonclassical properties and manipulation of the quan- 
tized electromagnetic (EM) field have been a center of 
interest in quantum optics jjj. A variety of methods 
have been proposed for the generation of squeezed states 
of the EM field, and several of them have been realized 
experimentally The quadrature-amplitude squeezed 
state in which fluctuations of the in-phase or out-of-phase 
component are suppressed to below those of the coher- 
ent state can be generated via nonlinear optical pro- 
cesses ||[|]. The photon- number squeezed state exhibit- 
ing sub-Poissonian photon statistics can be generated us- 
ing semiconductor lasers (5j, light-emitting diodes 
and tailor-made semiconductor heterostructures Q. 

A coherent state with average photon number n has 
the relative photon- number fluctuation of An/n — \j \ffi 
and the phase fluctuation of Aqb ~ l/(2^/n). Hence, if 
n S> 1, there is practically no need of squeezing. In 
a few-photon regime, however, An becomes comparable 
to h and A<f> becomes of the order of one. It is thus 
within the few-photon regime that the manipulation of 
quantum fluctuations becomes crucially important. Such 
quantum-controlled few-photon states might be useful, 
e.g., for optical interconnections in semiconductor mi- 
crostructures and spectroscopic diagnostics in biology. 

In the present paper, we develop a method of manip- 
ulating a collective atomic state to generate quantum- 
controlled few-photon states M. Radiation from atoms 



has been extensively studied in quantum optics |10|, e.g., 
superradiance [ pT| , resonance fluorescence |l2[ |, photon 
emission in the cavity | fi3| , etc. As regards nonclassical 
properties of radiation, it is known that resonant fluores- 
cence exhibits photon anti-bunching and sub-Poissonian 
photon statistics [^4 15 1. It was pointed out in Ref. |l6| ] 



that squeezing in the resonant fluorescence is related to 
quantum fluctuations in the atomic state. However, the 
relation between quantum fluctuations of the collective 
atomic state and those of the emitted photon state has 
yet to be fully explored from the standpoint of the con- 
trol of few-photon states. The aim of the present pa- 
per is to show that we can generate quantum-controlled 
few-photon states by preparing the atoms in a squeezed 
atom state (SAS), which is a collective state of quantum- 
mechanically correlated atoms whose quantum fluctu- 
ations in population difference or collective dipole are 
suppressed to below those of the coherent atom state 
(CAS) Q. The SAS can be generated via the inter- 
action of atoms with a coherent state of photons in the 
cavity having a high quality factor p8[ . It will be shown 
in Sec. fv| that the SAS can be used as a tunable source of 
squeezed radiation. This is owing to the fact that quan- 
tum fluctuations of the atomic state arc rather faithfully 
transferred to those of the emitted photon state. It will 
be shown that the number-phase uncertainty relation of 
photons can be manipulated only if the atoms are in the 
SAS. 

It is well known that the state of a two-level atom can 
be mapped onto that of a spin 1/2. A collection of N 
two-level atoms can be described with a system of spins 
whose magnitudes are at most N/2. In particular, if all 
atoms are in the same pure state, the collective atomic 
state can be described by a single spin N/2. The con- 
cept of squeezing in the spin or SU(2) algebra 
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provides a mathematical definition of squeezed states in 
a system of two-level atoms and in other systems that 
can be described by the spin algebra. Yurke et al. |2l]] 
has pointed out that the Mach-Zehnder interferometer 
is described by spin, and that its phase sensitivity can 
reach the fundamental limit of 2/N using an N particle 
squeezed state. Kitagawa and Ueda |25j showed that such 
a squeezed state can be realized using the Coulomb inter- 
action between charged particles. Wineland et al. |2(| ap- 
plied the SAS to Ramsey spectroscopy and showed that 
its sensitivity can surpass that of uncorrelated atoms. 

This paper is organized as follows. Section [fl] briefly 
reviews the interaction between two-level atoms and pho- 



tons in the cavity. Section HI defines the SAS in terms 
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of the spin representation of two-level atoms and dis- 
cusses its physical meaning. Section IV analyzes dy- 



namical processes to generate the SAS. Section [v| de- 
scribes how quantum-controlled radiation is generated 
from squeezed atoms. Section VI discusses possible ex- 



perimental schemes to implement our theory. Some com- 
plicated algebraic manipulations are relegated to the ap- 
pendices to avoid digressing from the main subjects. 



II. INTERACTION BETWEEN PHOTONS AND 
TWO-LEVEL ATOMS IN THE LOSSLESS 
CAVITY 

The EM-field operator in a lossless cavity can be writ- 
ten as 



E(r) 



[f„(r)a n -f*(r)4] , (1) 



where and d„ are the creation and annihilation op- 
erators of the EM field for the nth mode, and f,i(r) 
is the corresponding orthonormal mode function satis- 



fying (V 2 



\/c 2 )t n = 0, V • f„(r) = 0, and on the 



boundary the tangential component is required to van- 
ish: f„||(r) = 0. The Hamiltonian of the EM field in the 
cavity is given by 



(2) 



Here and henceforth, zero-point energies are ignored be- 
cause they do not affect the following discussions. 

Suppose that atoms have the upper energy band \ej a ) 
and the lower energy band \gjp), where j distinguishes 
atoms, and a and j3 denote Zeeman sublevels, if any, of 
the upper and lower energy bands, respectively. When 
the sublevels in each band are degenerate, the Hamilto- 
nian of N identical atoms have the form 



JY 



TlUJA 



\ e 3c}(e ja | - \9if>) (MP I I > ( 3 ) 




where Tiuja is the energy separation between the two 
bands. We consider a situation in which a collection of 
two-level atoms is placed in the cavity, and interacts with 
the EM field via the electric-dipole interaction described 
by 



N 



^/ = -^D 3 -E(R,), 



(4) 



where D, 



-e^ fc (fjfc — Rj) denotes the electric-dipole 
operator of the jth atom, which is the sum of differences 
between the position of the nucleus Rj and the posi- 
tions of the electrons f 3 i~ that belong to the jth atoms. 



We neglect the dynamics of the center-of-mass motion 
of atoms, and replace Rj with a c-number. Making the 
rotating- wave approximation in the Hamiltonian (|]) , we 
obtain 



N 



j = l a, ft n 



\^ n {Hj)\e ja )(e 3a \D 3 \g j p)(g jf} \a n - H.c. 



(5) 



where H.c. denotes the Hermite conjugate of the preced- 
ing term. 

We assume that only a single mode of the EM field 
having energy Jkop and a single state for each atomic 
energy band | ) and \gj) participate in the interaction, 
and omit the subscripts n, a, and j3 in the following dis- 
cussions. The Hamiltonian (||) then reduces to 

N 1 

H * = E 2 [ dj ' 5(Rj) &§j+ + d i ' f *( R i) fit V] . (6) 



where £ (Rj) = — iW2huj/sof (Rj) is the amplitude of 
the electric field per photon, d^ = (e 3 ;|Dj \gj) is the 
electric-dipole matrix element, and §j+ = \sj){gj\ and 
Sj- = \gj)(ej\ are the raising and lowering operators for 
the jth atom. Provided that the dipole moment is the 
same for all atoms, the subscript j in d^ may be omit- 
ted. We define three operators Sj X = (s J+ + Sj_)/2, 
Sjy = ~ sj-)/2i, and s jz = (\ej){ej\ - \gj){gj\)/% 
which can be verified to obey the spin commutation re- 
lation [sj X ,Sj'y] = i8jj'§j Z and its cyclic permutations. 
The two-level atom can therefore be described by spin 
1/2. The subscripts x, y, and z do not denote spa- 
tial directions, but the expectation value of the operator 
Sj z + i represents the probability of the jth atom being 
found in the excited state, and s 3X and §j y indicate the 
quadrature-phase components of the oscillating dipole. 
This can be seen by rewriting the dipole operator in the 
form 



t> J =d\e J )(g J \+d*\g J )(e J \ 



= ds 



2 [Re(d)s JrE — Im(d)s 



jyi 



(7) 



The spatial direction of the dipole depends on how we 
excite atoms. For example, if the electric field at the 
position of an atom is linearly polarized, the dipole oscil- 
lates along the same direction. If the electric field at the 
position of the atom is circularly polarized, the dipole 
also rotates in time. 

Suppose that all atoms are located in a region small 
in comparison with the wavelength of the field, but that 
they are not located too closely together in order to avoid 
direct interactions between them. The Hamiltonian of 
the entire system is then given by 
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H = hujpa) a + TlloaS z + fig \aS+ + S-j . (8) 

where the coupling constant g = £(Rj) ■ d/(2h) is taken 
to be real without loss of generality, and the collective 
spin operators are defined by 

^ = X] S iA< & = ^'f' 2 )' ( 9 ) 

3 

and S± = S x ± iS y . It is easy to show that these collec- 
tive operators follow the commutation relation of spin, 
[S x ,S y ] = iS z , and its cyclic permutations. The Hamil- 
tonian (||) is referred to as the Jaynes-Cummings (JC) 
Hamiltonian P7| . 

It is worth pointing out that one can introduce the 
collective spin operators when the magnitudes of £(Rj) ■ 
dj in the Hamiltonian (^) are the same for all the atoms 
but their phases are different due, e.g., to different spatial 
locations of the atoms. The collective spin operators in 
this case may be defined as 

S^^e^i, & = (10) 

j i 

where 4>j is the phase of £ (ftj) • dj. For example, when 
atoms are located in a one-dimensional standing wave at 
every half wavelength, we have <j>j — jir. The operators 
( |l0| ) also satisfy the spin commutation relations and the 
Hamiltonian of the system is given by Eq. (^) in which 
S± is replaced by S'± . Even if the spin state described by 
the operators ([Io|) and that described by the operators 
(0) are the same, the corresponding states of atoms are 
different. When atoms are located in the same place 
the dipoles oscillate in phase. When they are located at 
every half wavelength, the neighboring dipoles oscillate 
out of phase. Nevertheless, the photon states generated 
by these atoms via the JC Hamiltonian (J§|) are the same. 

When we move to the rotating frame for both the pho- 
ton field and the atoms via a unitary transformation 
U (t) = e i("F& f &+"ASz)t^ thc Hamiltonian (g) is trans- 
formed to 

BTJ 

UoHUl + »fi-0^#o = gfi(aS + e- M + a?S-e iSt ), (11) 

where 6 — ojf — ^>a denotes the detuning between the 
atoms and the field. When 5 is zero, Eq. (jll]) becomes 

H mt = hg (aS+ + a^S-^j . (12) 

This commutes with the rotation operator, 

U{ip) = e- iv ^' a+S '\ (13) 

and is therefore invariant under rotation. This rotational 
invariance allows us to choose a convenient frame of ref- 
erence without loss of generality. For instance, when ini- 
tially the EM field is in the coherent state \a) and the 



atoms are in the fully excited state \S,M — S), we can 
arbitrarily choose the phase of the initial coherent state 
without loss of generality. Time development from the 
other initial state \ae~ %v ) \S, M — S) can be obtained by 
a mere rotation U(ip). 

III. SQUEEZING IN COLLECTIVE TWO-LEVEL 
ATOMS 

As shown in the preceding section, a collection of two- 
level atoms can be described by collective spin operators 
<M). An eigenvalue of the Casimir operator S 2 = S% + 
+ Si is given by S(S + 1), where the total spin S 
can take on values, S = lj> y ~ 1, • • *, (or 1/2) when 
the number of atoms N is even (or odd). For each total 
spin S there are iV!(2S* + l)/[(|iV + S + l)!(|iV - 5)!] 
different subspaces. Generally speaking, a state of N 
two-level atoms can be described by a mixture of these 
subspaces. 

Because the JC Hamiltonian (§h is described by the 
collective spin operators which never mix the subspaces 
having different total spins, we will restrict our discus- 
sions to a single subspace having the maximal total spin 
N/2. This state can be most easily accessed from the 
state in which all the atoms are either in the ground state 
or in the excited state. It is interesting to note that the 
subspaces having the same total spin behave exactly the 
same within the JC model if the numbers of atoms are 
different. For example, the state of two atoms having the 
total spin 1 and that of 100 atoms having the same total 
spin 1 obey the same JC Hamiltonian. No single-mode 
photon field distinguishes between these atomic states 
through the JC interaction. 

A state of the single-mode photon field is defined as 
squeezed if, for a nonzero range of parameter <f>, ((Aa^) 2 ) 
is smaller than that of the coherent state — the standard 
quantum limit (SQL) — of 1/4, where is defined as 

a^ = ^{ae-^ + aU^). (14) 

The canonical commutation relation is given by 
[a0,d^ +7r y 2 ] = i/2, and the conventional in-phase and 
out-of-phase components a% and &2 can be expressed as 
Si = o.0 = o an d «2 = a0=7r/2i respectively. From the com- 
mutation relation we have 

((Aa,) 2 )((Aa 0+ . /2 ) 2 )>l. (15) 

The coherent state has the variance of ((Aa^) 2 ) =1/4 for 
any <f> and satisfies the equality in Eq. (|l5|). The profile 
of quantum fluctuations of a photon state described by 
a density operator pp can be visualized with the quasi- 
probability distribution 

Q(a) = -{a\p F \a), (16) 
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where \a) is the coherent state with amplitude a. The 
quasi-probability distribution of the coherent state is 
isotropic and that of the quadrature-amplitude squeezed 
state is elliptic. 

The coherent state of a spin-5 system is defined by 

\9, <j)) ee exp[i9(S x sine/)- S y cos^)]|5, M = 5) 
s 



E 

M=-S 



25 
S + M 

S-M 



a i(S-M)<f> 



x sm- 



S+M 



COS ■ 



\S,M), (17) 



which is referred to as the coherent spin state (CSS) 
or the Bloch state |Tij| . The mean spin vector of the 
CSS \6,(f>) points in the direction n = (S)/|(S)| = 
(sin cos 0, sin 9 sine/), cos 9), where |(S)| = ((S x ) 2 + 
(S y ) 2 + (Sz) 2 ) 1 / 2 . Denoting m as the unit vector that 
is normal to both n and the unit vector of the z direc- 
tion e z , namely, m = nxe z /|nxe z | = (sin cj), — cos </>, 0) , 
we may express \6,(f>) as 



\6,<t>) = exp[i0m-S]|5,M = 5). 



(18) 



When a system of two-level atoms is described by 
Eq. ([I?]) in the spin representation, we will say that the 
atoms are in a coherent atom state (CAS). The compo- 
nent of S normal to the mean spin vector is given by 

5(n, x) = exp(— i\S ■ n)(S • m) exp(z%S • n), (19) 

where \ denotes the angle defined on the plane normal to 
the mean spin vector (see Fig. 1(a)). The commutation 
relation between the two quadrature components is given 
by 

[S(n,x),S(n,x + 7r/2)]=iS-n, (20) 
and the corresponding uncertainty relation is given by 

l(S>| 2 



({AS(n,x)y){[AS(n, X + ir/2)y} > 



The CSS satisfies the equality in the uncertainty rela- 
tion @, and ([A5(n,x)] 2 ) = 5/2 for any X - The CSS 
therefore has an isotropic fluctuation normal to the mean 
spin vector as shown in Fig. ^(a), where the spin state is 
visualized with the quasi-probability distribution of spin 
defined by 



Qs 



25 + 1 
1 

47T 



■\Pa\ 



(22) 



where pa is the density operator of a collective atomic 
state. Analogous to the case of photons, a spin state is 
defined as squeezed if the following inequality holds for 
a certain x- 



S(n,5C) 




FIG. 1. The quasi-probability distributions of (a) the co- 
herent spin state and (b) the squeezed spin state. The unit 
vector n points in the direction of the mean spin vector, and 
the unit vector m is normal to both n and the S z direction. 
The spin component S(n, x) is normal to the mean spin vec- 
tor, and the angle x i s measured from m. 



<[A5(n, x )] 5 



< 



(23) 



That is, the squeezed spin state (SSS) is a state whose 
fluctuation of one component normal to the mean spin 
vector is less than half of the length of the mean spin 
vector. When the condition ( p3|) is satisfied, the variance 
(21) of the quadrature component ([A5(n, x + 7r /2)] 2 ) must 
be larger than |(S)|/2 in order to obey the uncertainty 
relation (pl|), and hence the fluctuation profile on the 
spin sphere becomes elliptic, as shown in Fig. 0(b). 

Squeezing in spin or angular momentum has been dis- 
cussed by many authors 0,|l^-^6ll. However, the def- 
initions of the SSS in Refs. 




depend on 



|6|,y,g0| 

the specific spin coordinates and are therefore not invari- 
ant under rotation in the spin space. It was pointed out 
in Ref. |25| that the direction of the mean spin vector 
n should be taken into account to define the SSS in a 
rotation-invariant manner as in Eq. (pa). 

Mathematically, SSS satisfying the condition ( |23| ) 
can be generated by unitary transformations from 
the CSS. The unitary transformations have the forms 
exp(-i£5j)|0 = tt/2,0) and expH^S^ -Sl)}\9 = 0,0), 
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where £ and r\ denote the parameters that characterize 
the degree of one-axis twisting and that of two-axis coun- 
tertwisting, respectively p5j . 

Let us return to the spin representation of two-level 
atoms. We define the squeezed atom state (SAS) as a 
state of two-level atoms that are in the SSS in the spin 
representation. We note that quantum-mechanical corre- 
lations between atoms must be established for the atoms 
to be in an SAS. The state in which all atoms are in their 
ground state is in a CAS \0 = ir, </>) in the spin represen- 
tation, and not in an SAS. If they are irradiated by a 
ir/2 pulse, the spin state becomes \8 = n/2,4>), which is 
also not squeezed, because atoms are described by the 
same state and are not quantum-mechanically correlated 
with each other. We also note that a single atom cannot 
be squeezed, since ([A5(n, x)] 2 ) is always l/4(= 5/2) 
for spin 1/2. In other words, the single atom cannot 
be squeezed because it has no partner with which to be 
quantum-mechanically correlated. 

According to the definitions of the collective spin oper- 
ators (||) , S z represents the population difference of two- 
level atoms, and S x and S y represent quadrature-phase 
components of the electric dipole. Squeezing of the S z 
component thus means reduced fluctuations in the pop- 
ulation difference at the expense of the enhanced dipole 
fluctuation, while squeezing of S x , S y , or their arbitrary 
linear combination 



^ = 2 (^ +e 



(24) 



means reduced dipole fluctuations at the expense of the 
enhanced fluctuations in the population difference. 

To measure the S z component, one can use an ioniza- 
tion detector which counts the number of atoms in the 
excited state. If such measurement is carried out repeat- 
edly, with the atoms prepared in the same state for every 
measurement, the variance of the population difference 
((AS Z ) 2 ) is obtained. Variances of the other spin com- 
ponents can be measured by rotating the spin state so 
that they become the S z component. The rotation in 
the spin space can be realized by irradiation of maser or 
laser with classical intensity to the atoms. The frequency 
of the maser or laser is assumed to be resonant with the 
transition frequency of the atom. The Hamiltonian de- 
scribing the irradiation process of the classical field is 
obtained by replacing the operator a with a c-number a 
in the JC Hamiltonian (Of), 



(25) 



Eel = hg(aS+ + a* S-) 

= 2Hg\a\(S x cos (f> c - S y sin0 c ), 



where <f> c = argq is the phase of the classical field. The 
Hamiltonian (25) rotates the spin vector about the axis 
through angle 2g|a|X'i, where Ti is the irradiation 
time. For example, the S x component can be measured 
by counting the population difference with the ionization 



detector after irradiation of the classical field correspond- 
ing to the operation exp(— i^S y ). In this operation the 
collective dipole of the atoms S x is converted to the pop- 
ulation difference S~. 



IV. PREPARATION OF SQUEEZED ATOM 
STATES 

Several schemes for generating the SAS have been pro- 
posed. Barnett and Dupertuis ||| considered the inter- 
action of the antisymmetric collective dipole with the 
coherent EM field. Agarwal and Puri p3 examined 
the steady state of atoms interacting with broadband 
squeezed radiation. Although a coordinate-dependent 
definition of spin squeezing ((AS x ( y )) 2 ) < |(5 z )|/2 is 
used in Refs. Q and Q, the states constructed there 
also satisfy the coordinate-independent condition (p2|). 
Wineland et al. |26| considered the stimulated Raman 
coupling between kinetic motion of atoms in an ion trap 
and internal levels of atoms, and showed that by ini- 
tially squeezing the kinetic motion one can generate the 
SAS of the internal state via the JC interaction. They 
also showed that the coherent state of the kinetic motion 
can generate the SAS via the parametric-type interac- 
tion. Kuzmich et al. p8j considered V-type three-level 
atoms driven by squeezed light that leads to the SAS. 

In the present paper we follow the scheme proposed 
in Ref. fl8|| , namely the interaction between the atoms 
and the coherent state of photons in a high-Q cavity. 
The higher-order interaction between atoms and photons 
establishes the quantum correlation between the atoms, 
thereby reducing the dipole fluctuation. This scheme is 
simple in that no special field state, other than the co- 
herent state, is required. 



A. Analysis for the case of two atoms 

The JC model can be solved exactly for up to three 
atoms, and in the zero-detuning case for up to eight 
atoms. We will henceforth assume zero detuning 5 = 0, 
and employ the Hamiltonian (^). By exactly solving 
the dynamical evolution for two atoms, we discuss the 
properties of this system. 

We consider the case in which both atoms are initially 
in the excited state \S — 1, M = 1) = \1,1)a and photons 
are in an arbitrary superposition state ^ n c n \n) F , where 
\n) p is the photon-number state. The time development 
is calculated to be [18| 

\ 1 p(t))=e-^ t ° H \n) F \l,l) A 

OO 

= Y / Cne' l{n+1) ^ t [pn(t)\l,l) A \n) F 
n=0 

+q n (t)\l,0) A \n + l) F + r„(t)|l, -l) A \n + 2) F ] , 
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(26) 



where 



Pn(t) = 



(n + 1) cos ^2(2/1 + 3)gt + n + 2 



2n + 3 



(27a) 
(27b) 



r n (t) = ^ {n ^ll + 2) (cos VWn + 3)gt - l) . (27c) 

One can calculate any physical quantities from this solu- 
tion. 

Let us first consider the photon- number state \n)p as 
the initial state. In this case the initial state |n)f|l, 1)a is 
invariant with respect to rotation (|l3|), and consequently 
(S x ) — (S y ) — 0, which remains true at later times. The 
variances of the components normal to the mean spin 
vector are calculated to be 




((A5,) 2 ) = <(A5,) 2 } 



yJ 

1 A , n + l . 2 
t: 1 + 7i ~ sm 

2 1 2n + 3 



FIG. 2. Time evolutions of the normalized variance 
2((A&) 2 ) and the squeezing factor 2{(A&) 2 )/| (S) | for two 
atoms. The two atoms are initially excited, and the EM field 
is initially in the coherent state with amplitude a = 10. The 

i— \ solid curves show the numerical results, and the dashed ones 

V2(2n + 3)gt j , (28) ghow appi . oximate solutions @ and @. 



which is always greater than 5/2 = 1/2, and hence 
the spin state can never be squeezed. Generally, when 
the initial state is invariant with respect to the rotation 
U(ip), the atoms can never be squeezed for any number 
of atoms. 

When the photon field is initially in the coherent state 
|a), the coefficients are given by c„ = e' a ' < 2 a n /vri\. 
The amplitude a can be taken to be real without loss 
of generality, and in this case (02) and (S x ) vanish at 
any time (see appendix . Therefore the S x direction is 
always normal to the mean spin vector. The variance of 
S x is calculated to be 



+i (sm2^/ngt + sin 4-\/ngt) , 
8n \ / 



Sn 

.-(9*) 2 /2 



(31a) 



cos 



2\fngt 



5gt_ 



Sill 



2\fngt 



sin 2 1\fnqt. 



(31b) 



Therefore, if the squeezing factor defined by 

((AS^) 2 ) ( 9 *) 2 1 . r- 3 . ,„ r- 

- ; ' ' ~ e~^~ sin 2 y/ngt H sin 2 2>/figt 

I (S) |/2 h y 8n 

, 3 5 i . 



v2\fngt 



(32) 



<(A5 X ) 5 



- P -W\ 



CXJ 

y a 



2n— 2 



=p n (t)r n - 2 (t) 



+^"W 2 £ £L_ [p„(<) 2 + 2g„(t) 2 + r n (tf] 



n=0 



(29) 



When a» 1, the photon- number distribution has a nar- 
row peak relative to the mean photon number n, and one 
can expand ( p9| ) with respect to n — n. Replacing the 
summations with the integrals we obtain, for gt < 1/ 1/^, 



((A5,) 2 



— sin 4 + -^= sin 2\/%t. (30) 
2n 2\/n 



Similarly, (S y ) and (5 Z ) are approximated to be 

(S y ) ~ -e'(f*) 2 / 2 sin 2 v^fft - 4= (7 - ~ sin 2 v^i 

yn \4 2 



is less than one, the condition for the SAS (|23j) is ful- 
filled. Figure 2 compares the time evolution of the ap- 
proximate formula ( [32] ) (dashed curve) with the exact one 
which is numerically calculated from ( p6|) (solid curve) 
for two atoms and for n = a 2 — 100. The parame- 
ter gt in Fig. 2 and all the quantities appearing in the 
figures presented henceforth are dimensionless. We find 
that both curves are in excellent agreement and the SAS 
is attained around gt = 0.2. The variance of another 
component that is normal to both the mean spin vector 
and the S x direction never reduces to below 1/2. It can 
be shown numerically that the SAS never occurs after 
the first minimum around gt = 0.2. Although in Fig. 2 
the second minimum of the variance ((A5;,;) 2 ) goes be- 
low the first minimum, the squeezing factor does not go 
below the first minimum because the length of the spin 
vector also decreases. 



Equation (32) shows that squeezing vanishes when the 
intensity of the coherent state is sufficiently large, fj>l, 



6 



which is due to the fact that the classical field merely 
rotates the spin vector. The photon-number state cannot 
produce the SAS, as mentioned above. We thus find that 
both wave and particle aspects of photons are necessary 
for atoms to be squeezed. 



B. Analytic approach for the case of a large number 
of atoms 

We provide here approximate analytic expressions for 
the case of a large number of atoms. These are derived 
by neglecting the terms of order 1/N relative to the dom- 
inant terms in the equations of motion, which are there- 
fore very accurate when the number of atoms N is very 
large. 

The initial state is assumed to be the totally excited 
state of the atoms \S,M = S) and the coherent state of 
the field |a), where a is assumed to be real and hence 
(S x ) = (0.2) = 0. The other averages obey the equations 
of motion (see appendix M for derivations), 



d(Sy) 

dt 

d(S z ) 

dt 
d(di) 

dt 



-2.g(a 1 )(^ 2 ), 



(33a) 
(33b) 
(33c) 



which become those of a pendulum, if we set 



(Sy) 


N ■ a 

= — suit/, 

2 


(34a) 


<&> 


N a 
= — cos 0, 

2 


(34b) 


(Si) 


1 d6 
~ 2g~dl' 


(34c) 



The solutions of Eqs. (33a) can be expressed in terms of 
Jacobi's elliptic functions (29). Solving the equations of 
motion for fluctuations, we obtain 



<(Aa 2 ) 2 ) 
((A^) 2 ) 



1 



4dn (u\m) 
,2, 



[l + mE 2 (u\m)] , (35a) 



N I sn (u\m)cn 2 (u\m) 
— < m 



+dn(zt|m) 



dn (u\m) 
sn(w|m)cn(u|m) 
dn 2 (it I m) 
2 ] 



E(u\m) 



(35b) 



where u = gt\/N + a 2 and m = N/(N + a 2 ). Jacobi's 
elliptic functions [g9[ are defined by sn(w|m) = sinp, 
cn(u|m) = cos tp, dn(u|m) = \fl — m sin 2 ip, where u and 
tp are related by 
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gt 

FIG. 3. Time evolutions of (S z )/S, 2((AS x f)/S and 
((Aa2) 2 ) for 100 atoms (S = 50). All the atoms are initially 
excited, and the EM field is initially in the coherent state 
with amplitude a = 10. The solid curves show the numerical 
results, and the dashed ones show approximate solutions j3^ ) 
and (B^a). 



de 



yl — rn sin 2 9 



(36) 



The elliptic integral of the second kind is given by 
E{u\m) — J Q " dn 2 (u'|m)du'. 

Figure ^compares the analytic solutions (|35| ) (dashed 
curves) with the numerically exact ones (solid curves) for 
100 atoms and a = 10. We find that the analytic curves 
are in excellent agreement with the numerical ones. The 
analytic curves, however, begin to deviate from the nu- 
merical ones at around gt ~ 0.3. This is because the dif- 



ferential equations (B4) and (B9) include errors of order 
1/N relative to the dominant terms, which accumulate to 
produce errors in the solutions of order e gt ^ /N, which 
becomes of order unity around gt ~ 0.3. 

The analogy to the pendulum gives us a quali- 
tative and simple account of the squeezing mecha- 
nism. When the pendulum points in the direction 
(sin 9 cos 4>, sin 8 sin <j>, cos 9) , it undergoes the force to- 
ward the direction (cos 9 cos </>, cos 9 sin 0, — sin 9). In the 
present case, where a is taken to be real, the pendulum 
begins to fall toward the negative S y axis and rotates 
Suppose that the pendulum has 
= -f + 5<f)), the 
direction of the force is (cos 08<fr, 0, — sin 9). This force 
increases the deviation when cos 9 > 0, and decreases it 
when cos# < 0. In fact, in Fig. |[ {^(AS^) 2 ) increases 
when (S z ) > 0, and decreases when (S z ) < 0. 



on the Sy-S z plane 
a deviation from the S y -S z plane 



7 



0.5 1 1.5 

gt 



20 40 60 

gt 



FIG. 4. Time evolutions of (S z )/S, 2((AS X ) 2 ) /\(S)\, (a) 
and {(An) 2 ) /(h) for 10 atoms (S = 5). All the atoms are 
initially excited, and the EM field is initially in the coherent 
state with amplitude a = 3.3. 



FIG. 5. Long-term behaviors of (S z )/S, ((AS X ) 2 ) and 
((Ah) 2 ) / (h) for 10 atoms. All the atoms are initially excited, 
and the photon field is initially in the coherent state with 
amplitude a — 3.3. 



C. Numerical analysis 

When the number of atoms is intermediate, analytic 
solutions are unavailable, so we study the dynamical evo- 
lution of the system by numerically diagonalizing the 
Hamiltonian (|l2|). The amount of computation increases 
with increasing the number of atoms N roughly as TV 3 . 
The initial state is assumed to be the totally excited state 
of the atoms \S,M = S) and the coherent state of the 
photon field \a), where a is again taken to be real. 

Figure [| shows time evolutions of statistical proper- 
ties of atoms and photons. The number of atoms is 10, 
and the amplitude of the initial coherent state is cho- 
sen to be a = 3.3 to obtain the maximal squeezing of 
the atoms. The S x component is always normal to the 
mean spin vector, since (S x ) = 0. In Fig. ||, the squeez- 
ing factor 2((AS X ) 2 ) /\ (S) | becomes less than one, which 
indicates that the SAS is obtained. The maximum de- 
gree of squeezing is attained in the first minimum. It 
is found from the long-term behavior that the squeezing 
never occurs at a later time. The fluctuation of the other 
component that is normal to both the S x direction and 
the mean spin vector never fulfills the squeezing condi- 
tion (p3[). Since the mean spin vector rotates in the S y - 
S z plane, (S z ) oscillates with the amplitude of | (S) |. The 
amplitude of the photon field also oscillates with the same 
period but out of phase because of the energy exchange 
between the atoms and the photon field. The variance 
((AS X ) 2 ) increases when (S z ) > 0, and decreases when 
(S z ) < 0, as discussed in the previous subsection. The 
Fano factor ((An) 2 ) /(h) of the photon field also goes be- 
low the SQL, and its behavior is very similar to that 



of ((AS X ) 2 ). The long-term behavior of this system is 
shown in Fig. [5]. The collapse and revival phenomena oc- 
cur in the population difference and in the Fano factor as 
in the case of a single atom |3(J . The revival peak of the 
Fano factor splits and there is a small revival before the 
main revival. The variance ((AS X ) 2 ), on the other hand, 
oscillates with the same period as the revivals, and the 
initially regular oscillations gradually change to random 
fluctuations around some value. 

The degree of squeezing of the SAS depends on the 
number of atoms N, and for each N the maximum de- 
gree of squeezing is attained at a particular amplitude a 
of the initial coherent state. Figure || shows the mini- 
mum squeezing factor 2((AS X ) 2 ) /\(S)\ for each number 
of atoms and the amplitude of the initial coherent state 
that gives this factor. We find that the higher degree of 
squeezing can be obtained for the larger number of atoms. 
The squeezing factor tends to behave as N~ - 25 when N 
is more than about ten, and the optimal amplitude a 
behaves as A^ 0,29 . 

V. QUANTUM-CONTROLLED RADIATION 
FROM SQUEEZED ATOMS 

It is natural to expect that the atoms whose collective 
dipole or population difference is squeezed can radiate 
the photon field having nonclassical properties. We will 
show that this is indeed the case, and that quantum fluc- 
tuations of the photon field can be controlled by manip- 
ulating the SAS, which is done by applying a classical 
field to the atoms. 
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FIG. 6. Minimum values of the squeezing factors 
2((ASx) 2 ) /\ (S) | obtained by the interaction of atoms with 
the coherent states of photons as a function of the number 
of atoms N. For each N the amplitude of the coherent state 
a is chosen to give the best squeezing factor. The squeezing 
factor tends to scale as jV -0,25 for large N and the optimal 
amplitude a as A r °' 29 . 
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Atom Field 

FIG. 7. The relation between (a) the quasi-probability 
distribution of the prepared atoms (G2) and (b) that of the 
emitted photon field (|l6|) . The angle <f) (or —(f> — tt/2) repre- 
sents the direction of the mean amplitude of the photon field 
(or the mean spin vector), and ip represents the direction of 
the fluctuations of the spin and the photon field. The S z 
component of the mean spin vector is negative. 



photon field are not correlated. Since the first derivative 
vanishes at t = 0, the time development for small t is 
determined by the second derivative. From Eqs. j3^a) 
and (H^b) we have 



A. Radiation from squeezed atoms 

The Heisenberg equations of motion for 
>S-0-7r/2 are written as 



[H lot ,a^] — gS- ( f,- 7r /2, 

2ga <t> S z 



rot q 

<3-<£-7r/2 



, and 

(37a) 
(37b) 



Equation (|37|a) indicates that the phase of the photon 
field is connected with the direction of the spin vector. 
When the spin vector is tilted toward the direction of 
—(f) — 7r/2, the field is initially amplified toward the di- 
rection of <p, as illustrated in Fig. |?]. The equations of 
motion for various fluctuations are given by 



dt 



((Aa^)(AS_^_ n/2 )) 



-((AS-^/j,) 2 ) 



2 ff ((Aa^)(A5 )), (38a) 

g ((AS_^ /2 ) 2 ) 

+2{(Aa i ,)(Aa i ,S z )) , (38b) 



2g([AS_ il _ 7T/2 ,Aa i ,S z }+), (38c) 



where AO = O — (O), and [A,B] + = AB + BA is an 
anti-commutator. The angle ip in Eqs. ( |38| ) represents the 
direction of the fluctuations of the spin and the photon 
field, as shown in Fig. The right-hand side of Eq. (|3"§|a) 
vanishes at t = 0, because initially the atoms and the 



^((Aa^) 2 ) = 2 5 2 [((AS_^ /2 ) 



2((Aa 4 ,)(Aa i ,S z )) 
(39) 



At t = 0, the right-hand side of Eq. (|39| ) reduces to 
2g 2 [((AS^ 7r/2 ) 2 ) + (S z }/2] because ((Aa^) 2 ) = 1/4 
for the vacuum state. Therefore, if the initial spin state 
satisfies the condition, 



{{AS^ /2 f) < 



(40) 



the photon field will evolve into a squeezed state. To 
satisfy the inequality (HQ ), (S z ) must be negative. The 
equation of motion (|39j) indicates that the fluctuation 
profile of the photon field is connected with that of the 
spin state. From Eqs. j3^a) and (|39|), then, the direction 
toward which the spin vector tilts corresponds to the di- 
rection of the displacement on the complex-a plane of the 
photon field, and the squeezed or enhanced direction of 
the spin fluctuation corresponds to that of the fluctuation 
of the photon field. Consequently, the quasi-probability 
distribution of the photon field on the complex-a plane is 
expected to behave like the quasi-probability distribution 
of the atoms on the spin sphere, as illustrated in Fig. [?]. 

When the tilting angle of the spin vector from the z 
axis is small, i.e., 9 ~ tt, we can approximately solve the 
equations of motion ( p7| ) and (|3^). In this case, (S z ) is 
almost constant, and S z can be replaced by a constant c- 
number (S z )q, where (• • -)q denotes the expectation value 
with respect to the initial state. With this approxima- 
tion, Eqs. (p7|) can be solved, giving 
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= (5 ° ^" sin-, _' - ^ i-llal 
2|(^)ol 



{S-4,-*/a) = (S-4-*/2) coBy/2\{S z )o\gt. (41b) 

The equations of motion for the fluctuations J38| ) become 
closed forms in this approximation, and the solutions are 
given by 



((Aa^) 2 ) = icos 2 y/2\(S 2 ) \gt 

— 1p — lt J 2) 2 )o . 



sin 2 \/2|(&) | 3 t, 
(42a) 



2|(^)o| 



{{AS-^/aY) = ((A5_^ /2 ) 2 ) cos 2 V /2|(&)o|.gt 

K^)ol 



+i^^ sin 2 p\(S z ) \gt. (42b) 

We find that if the condition ( ]4C| ) for the initial spin 
state is fulfilled, the variance of the quadrature am- 
plitude ( [42ft) goe s below the SQL of 1/4. At time 

t = it{2\J 2|(S' z )o|.g)~ 1 , the fluctuation ((Aa^,) 2 ) attains 
its first minimum 



((A«v) 



2 _ ((A5_^,_ 7r / 2 ) 2 )o 



(43) 



2|(^>o| 

and at the same time the amplitude of the field becomes 



maximum 



2|(^)o| 



(44) 



Figure g shows time evolutions of the amplitude and 
the variance of the photon field, where the initial atomic 
state is the SAS of 100 atoms. This atomic state 



is prepared by the method discussed in Sec. IV (the 
state at gt — 0.14 in Fig. ||). Since the tilting angle 
tan _1 (— (Sy)o/ (S z )o) — 0.258 is small, the small-angle 
approximation is valid. The solutions ( f4l| ) and (E3) are 
used for the theoretical curves in Fig. || (dashed curves). 
One can see that the analytic results agree well with the 
numerical ones (solid curves), and ((Aa 2 ) 2 ) goes below 
the SQL of 1/4. It can be shown numerically that the 
second and the later minimums of ((Aa 2 ) 2 ) are larger 
than the first minimum, and hence we should switch off 
the interaction when the first minimum is reached. 



B. Tailor-made radiation from squeezed atoms 

As illustrated in Fig. ^, the quasi-probability distribu- 
tion of the emitted photon state is like a projection from 
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FIG. 8. Time evolutions of the normalized variance 
2((ASx) 2 ) I S of the atoms, and the amplitude (a) and the 
variance {(A.CL2) 2 ) of the photon field emitted from them. The 
atomic state at gt = 0.14 in Fig. ^| is used as the initial atom 
state. The field is initially in the vacuum state. The solid 
curves show the numerical solutions and the dashed curves 
show the approximate solutions (fi"l|a) and (fl^). 

that of the prepared atomic state. This observation, to- 
gether with the solutions (41) and (pq), suggests to us 
that we can manipulate the direction of displacement and 
the direction of squeezing of photons by controlling the 
spin vector of the SAS. The rotation of the spin vector 
about an axis on the S x -S y plane can be made by apply- 
ing maser or laser with classical intensity to the atoms as 
described by the Hamiltonian (p5|). The rotation about 
the S z axis is realized by applying a dc magnetic field 
which causes a temporal detuning by the Zeeman shift. 
Combining these two processes, we can manipulate both 
the spin vector and the direction of squeezing. By ma- 
nipulating the SAS in the spin space, we can control the 
uncertainty ellipse of the photon field on the complex-a 
plane. Figure |^ shows the quasi-probability distributions 
of 100 atoms (left panels) and those of the emitted pho- 
ton states (right panels). In Fig. ^(a) the CAS is used, 
and in Figs. ||(b)-(d) the atom states are prepared in the 
SASs by the method discussed in Sec. IV, where the pa- 
rameters are optimized to obtain the maximum degree 
of spin squeezing (a = 6.8, gt = 0.19). The tilting angle 
of the spin vector from the negative S z axis is taken to 
be 7r/4 in Figs. ||(a)-(d), and the uncertainty ellipses are 
turned around by 0, 7r/4, and 7r/2 in Figs. ||(b), (c), and 
(d), respectively. One finds that the fluctuation profiles 
of the atomic states are rather faithfully transferred to 
those of the emitted photon states. Figures |(c) and (d) 
suggest that not only amplitudes and fluctuations but 
also higher-order moments of atom states are transferred 
to those of the photon states. We have thus demonstrated 



10 




8 



go 



-8 




8 8 

Re a 





FIG. 9. The quasi-probability distributions of 100 atoms 
(left) and those of the photons emitted from the atoms (right). 
In (a) the atoms are prepared in a coherent atomic state, and 
in (b), (c), and (d), they are prepared in squeezed atom states. 
In (b), (c), and (d) the uncertainty ellipses are turned around 
by angles 0, 7r/4, 7t/2, respectively. The mean spin vectors 
are tilted by 7r/4 from the negative S z axis. The spin sphere 
is seen from the negative S z axis. 



that by manipulating the SAS, we can control quantum 
statistical properties of the photon field at our disposal, 
which we would like to call tailor-made radiation. 

The squeezing of photons in the direction of phase can 
be obtained only if the atomic state is squeezed in the 
azimuth direction as in Fig. ||(b). Although the CAS 
can produce the photon-number squeezed state |nj as 
in Fig. ||(a), where the Fano factor is 0.81, it never pro- 
duces the phase-squeezed photon state by any rotation 
on the spin sphere. This can be verified numerically, and 
can also be deduced from the fact that the projection 
of the fluctuation profile on the complex-a plane from 



the spin sphere can never be squeezed in the direction of 
the phase if the fluctuation profile on the spin sphere is 
isotropic. To produce not only the amplitude-squeezed 
state but also the phase-squeezed state, the atom state 
must therefore be squeezed in the sense of the definition 



C. Available range of the tailor-made radiation 

Let us discuss the range of photon squeezing that is 
available by our method. We use the SAS generated by 
the interaction between the totally excited atoms and 
the coherent state of the photon field with an optimum 
amplitude as discussed in Sec. IV. The available range 
of the emitted photon field is obtained by plotting time 
evolutions of the radiation processes for various initial 
tilting angles of the spin vector of the SAS. 

Figure |l(] shows time evolutions of the amplitudes and 
the variances of the quadrature amplitudes of the pho- 
ton states emitted from the SASs of 100 atoms. Each 
trajectory is drawn with the initial tilting angle of the 
mean spin vector at every 7r/30. In Fig. |l0|(a), the SASs 
are prepared in the states squeezed in the longitudinal 
direction, as in Fig. ^(b) . The emitted photon states are 
therefore out-of-phase squeezed states. In Fig. |l(](b), the 
initial SASs are squeezed in the latitudinal direction as 
in Fig. |9|(d), and the emitted photon states are there- 
fore in-phase squeezed states. We find that in Fig. |lC|(a) 
the trajectories tend to return the same paths, whereas in 
Fig. [l0](b) the trajectories tend to round downward. This 
indicates that in the case of in-phase squeezing the en- 
ergy exchange and the fluctuation exchange between the 
atoms and the photon field tend to occur synchronously, 
and in the out-of-phase squeezing the fluctuation ex- 
change tends to be delayed against the energy exchange. 
When we draw the overlap region of Figs. [l(](a) and 
(b), we can obtain the available range of the quadrature- 
amplitude squeezed state. It can be shown that the larger 
number of atoms can produce the wider range of | (a) | and 
((Aa,/,) 2 ) ||. This is due to the fact that the larger is the 
number of atoms the larger will be the degree of squeez- 
ing of the SAS, as shown in Fig. ^. 

The ranges of the average photon number (h) and the 
Fano factor ((Ah) 2 ) / (n) available from the SASs and the 
CASs of 50 and 100 atoms are shown in Fig. O. It is 
found that for a given number of atoms the SAS can 
suppress the photon-number fluctuation more effectively 
than the CAS. The range of 100 atoms does not cover 
that of 50 atoms in Fig. [H} The SAS of 50 atoms can 
produce photon states having smaller Fano factors than 
the SAS of 100 atoms when the average photon number is 
less than about 40. For a given average photon number, 
therefore, there is an optimal number of atoms to produce 
the best photon-number squeezed state. 
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FIG. 10. Time evolution of the amplitudes and the vari- 
ances of the quadrature components of the photon state emit- 
ted from the squeezed 100 atoms. Each trajectory is drawn 
with the initial tilting angle of the mean spin vector at ev- 
ery 7r/30. The squeezed atom states are prepared in the 
same manner as in Fig. |^ and rotated to the states which are 
squeezed in the latitudinal direction in (a) (as in Fig. ^|(d)), 
and in the longitudinal direction in (b) (as in Fig. ^(b)). The 
dashed curves delimit the regions that the trajectories can 
reach. 



The ranges of the average photon number (n) and the 
phase fluctuation ((A<fi) 2 ) available from the SASs of 50 
and 100 atoms are shown in Fig. [l2]. Here we use the 
phase operator proposed by Pegg and Barnett |32| . When 
(a) is real and positive, the variance of the phase is ex- 
pressed as 



<(^) 2 > = T + E 



2(-l)™-" 
(n — n') 2 



F {n \pF\n)F, 



(45) 



where pp is the density operator of the photon field and 
\n)F is the photon- number state. Figure [l^ shows the 
variance of the phase ((A</>) 2 ) normalized by that of the 
coherent state having the same average photon number 
((A0) 2 ) co h. Here the phase is defined as squeezed when 
((A0) 2 )/((A</>) 2 ) coh is below unity. The dotted curve in 
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FIG. 11. Ranges of the average photon number (n) and 
the Fano factor ((An) ) / (n) of the photon field that can be 
obtained by the squeezed atom state (SAS) and the coherent 
atom state (CAS) of 50 and 100 atoms. The SASs are pre- 
pared in the same manner as in Fig. ^[ The regions above 
the curves show the available photon states. The solid curves 
show the results of the SASs and the dashed ones show those 
of the CASs. 
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FIG. 12. Available ranges of the average photon number 
(h) and the phase fluctuation ((A0) 2 ) normalized by that of 
the coherent state having the same average photon number 
((A^) 2 ) co h. The regions above the curves show the photon 
states that can be obtained by our method. The squeezed 
atom states of 50 and 100 atoms are prepared in the same 
manner as in Fig. K)l The dotted curve shows the lower bound 
of ((A0) 2 >/((A0) 2 > coh of the photon field. 



Fig. |l2| shows minimum values of ((A0) 2 )/((A(^) 2 ) co h for 
given average photon numbers, which are obtained by the 
method of Lagrange multipliers |35| l (see appendix |^). 
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The range of 100 atoms does not completely include that 
of 50 atoms as in the case of the Fano factor, which indi- 
cates that for a given average photon number there is an 
optimal number of atoms to reduce the phase fluctuation. 

In experimental situations, loss of photons in the cav- 
ity and spontaneous emission of atoms are unavoidable, 
and we therefore evaluate how much cavity loss and spon- 
taneous emission are allowed in order not to destroy the 
squeezing of the atoms and that of the photon field. We 
adopt the master-equation approach to take into account 
the effects of dissipation. The master equation in the 
presence of cavity loss and spontaneous emission is given 

by[@ 



dp 

at 



If 



(2apa) — a) ap 



h Lrl " ' ' 2 

(^2S-pS+ — S+S-p — pS+S- 



pa^a) 



(46) 



where p denotes the density operator of both the atoms 
and the photon field, and jj 1 and 7" 1 are the lifetimes 
of a single photon and a single atom in the cavity. We 
obtain time evolution of the density operator by numeri- 
cally integrating the master equation ( [i"6|) by the Runge- 
Kutta method. Figure |l3|(a) shows the contour plot of 
the minimum attainable values of the squeezing factor 
2{(AS X ) 2 ) /\(S)\ of the SASs obtained by the interaction 
of 10 atoms with the coherent state of the photon field. 
The amplitude of the coherent state is optimized to ob- 
tain the maximum degree of squeezing for each jf and 
7 a . Figure |l3|(b) shows the contour plot of the mini- 
mum values of ((Aa^) 2 ) of the photon field emitted from 
the squeezed atoms prepared in Fig. |l3|(a). The param- 
eters 7/ and 7 a in the radiation process are assumed to 
have the same values as in the preparation of the SAS. 
These results show that the generation of the SAS and 
the squeezed radiation are possible even in the presence 
of dissipation in experimentally feasible situations. We 
will discuss some concrete numbers in the next section. 



VI. POSSIBLE EXPERIMENTAL SITUATIONS 

We discuss possible experimental situations to imple- 
ment our theory. Our procedure of generating quantum- 
controlled few-photon states consists of three stages: (1) 
preparation of the SAS, (2) manipulation of the SAS (ro- 
tation of the spin vector in the spin space), and (3) radi- 
ation from these atoms. 

A simplest realization of our theory would be to fly a 
bunch of atoms through two cavities and a waveguide as 
schematically illustrated in Fig. |lj. This type of experi- 
ment may be done in a microwave regime, since the atoms 
are required to be within a region much smaller than the 
wavelength. If we use, for example, the 63p 3 / 2 <-> 61d 3 / 2 
transition of rubidium atoms, the resonant frequency is 
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FIG. 13. (a) The contour plot of the minimum values 
of the squeezing factor 2{{AS X ) 2 )/\{S)\ of the SASs obtained 
by the interaction of 10 atoms with the coherent state of the 
photon field. The amplitude of the coherent state is optimized 
to obtain the maximum degree of squeezing for each 7/ and 
7a. (b) The contour plot of the minimum values of ((Aa^>) } 
of the photon field emitted from the squeezed atoms prepared 
in (a). 



21.5 GHz, the wavelength is A ~ 10~ 2 m, and the cou- 
pling constant is g ~ 10 4 Hz. First, an atomic beam from 
an oven is collimated and velocity-selected. The variance 
of the velocity of the atoms must be Av <C A/T ~ 10 2 
m/s, where T is the time it takes the atoms to pass 
through the apparatus. A mechanical shutter can pre- 
pare a bunch of atoms from the atomic beam. The atoms 
in the bunch are then excited to the Rydberg state that 
is the upper state of the relevant two-levels, and enter 
the first cavity in which the photon field is prepared to 
be in a coherent state \a). The SAS is generated there 
by the higher-order interaction of the atoms with the co- 
herent state. The interaction time is gti ~ 10" 1 , i.e., 
ti ~ 10~ 5 s, e.g., in the situation in Fig. [|. The veloc- 
ity of the atoms is therefore required to be v ~ 10 3 m/s. 
The atoms then pass through a waveguide, where the 
atoms are irradiated by a pulse of microwave with clas- 
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FIG. 14. Schematic illustration of an experimental setup 
to implement the tailor-made radiation. The state of the 
atoms at each stage is shown with the spin quasi-probability 
distribution. A bunch of two-level excited atoms that is in a 
coherent atom state (CAS) is prepared by an oven, a colli- 
mator, a velocity selector, a shutter, and a pulse that excites 
the atoms. The atoms then go into the first cavity and inter- 
act with a coherent state of the photon field \a). The output 
atoms are in a squeezed atom state (SAS) . By the interaction 
with a microwave pulse in a waveguide, the mean spin vector 
is rotated to a desired direction. The atoms finally go into 
the second cavity and emit photons there. Left in the second 
cavity is the desired few-photon state. 



sical intensity, by which the spin vector representing the 
state of collective atoms is rotated. To control the rota- 
tion axis of the spin vector, the relative phase between 
the microwave and the coherent state in the first cavity 
must be controlled. The irradiation time of the classical 
field is much shorter than the interaction time in both 
cavities. Finally, the atoms pass through the second 
cavity in which the desired state of photons is emitted 
from the atoms. The interaction time is gt 2 ~ 10 _1 , i.e., 
i 2 ~ 10~ 5 s, e.g., in the situation in Fig. ||. The atoms 
thus pass through the two cavities within a few periods of 
duration 10 _5 s, which is much shorter than the lifetime 
of the Rydberg atoms ~ 10~ 3 s and the cavity lifetime 
~ 10 -1 s p5fl . From Fig. [ll| this cavity lifetime corre- 
sponding to jf/g ~ 10 -3 does not affect the squeezing. 
If the circular Rydberg states are used, the lifetime is ~ 1 
s, and decays from the relevant levels become negligible. 
Since the microwave frequency is used, the temperature 
should be lower than ~ 1 K in order to make the average 
number of thermal photons in the cavity much smaller 
than that of the produced photons. 



Another possible scheme is to use atoms confined in 
an ion trap or a magnetic trap in which the quantized 
kinetic motion of the atoms replaces the role of photons. 
Wineland et al. [^6| proposed the JC interaction between 
the Zeeman doublet of electronic states of each ion and 
the center-of-mass (CM) motion of an ensemble of ions 
via the inhomogcneous magnetic field. They pointed out 
that the stimulated Raman transition can also be used to 
couple the internal states of each ion to the CM motion 
of ions p6| , p6| . In these models the operators a and a) 
in the JC Hamiltonian (JsJ) are not for photons but for 
the quantized CM motion of ions in a harmonic trap. By 
using the stimulated Raman technique, our theory might 
be implemented as follows. First, the internal levels of 
the trapped ions are excited and the CM motion is cooled 
to the ground state ^] . The CM motions of two orthog- 
onal directions, say the z and x directions, correspond to 
the photon fields in the first and the second cavities in 
the method discussed in the preceding paragraph. In the 
first stage, the coherent state of the CM motion in the z 
direction is prepared and the Raman beams in this direc- 
tion are applied. The coherent state of the CM motion 
can be generated by sudden displacement of the trap cen- 
ter. When the atomic internal state becomes the SAS, 
the Raman beams are switched off. In the second stage, 
the Raman beams that do not affect the CM motion are 
applied, which rotate the spin vector in the spin space. In 
the third stage, the Raman beams in the x direction are 
applied, and the internal states of the ions are coupled 
to the CM motion in the x direction. By this coupling 
the information of the internal states is transferred to 
the CM motion in the x direction, which may be called 
a tailor-made motional state. Although this is not radia- 
tion, the method using the trapped atoms might be used 
to test our theory. 

The use of dielectric spheres as optical cavities might 
be another possibility, where the optical whispering 
gallery (WG) mode in the microsphere is employed. With 
the microsphere cavity, very low threshold lasing has 
been observed J38|,|3?|, and the Q value of more than 10 9 
has been achieved with highly transparent silica glass 
optical- fiber material jfO|. The atoms are fixed on the 
substrate and are coupled to evanescent waves of two 
microspheres which are placed very closely. They have 
the slightly different resonant frequencies uj\ and u>2- The 
optical WG mode in the first microsphere is prepared in 
a coherent state, while that in the second one is pre- 
pared in the vacuum state. In the first stage, the atoms 
are brought into resonance with the WG mode frequency 
in the first microsphere u)\, and are far from resonant 
with that of the second microsphere lo 2 - This can be 
done by Zeeman-shifting the transition frequency of the 
atoms by a magnetic field. When the atomic state be- 
comes the SAS by interacting with the coherent state, the 
interaction with the WG mode in the first microsphere 
is switched off by switching off the magnetic field. In 



14 



the second stage, the spin vector is rotated by applying 
a laser pulse resonant with the transition frequency of 
the atoms. In the third stage, by applying an appropri- 
ate magnetic field, the atoms are brought into resonance 
with the WG mode of the second microsphere u>2- By 
switching off the magnetic field, the desired photon state 
is left in the second microsphere. The coupling constant 
g can be of order 10 s , and 7//g ~ 10 -2 and "f a /g ~ 10 -1 , 
where the spontaneous emission rate of an atom in the 
free space is assumed. From Fig. [D| we find that both 
the SAS and the squeezed photon state are not washed 
out by the effects of dissipation. 



VII. CONCLUSIONS 

In conclusion, we have shown that quantum fluctua- 
tions of few-photon states can be controlled by using the 
SAS. This controllability is based on the fact that quan- 
tum fluctuations of the atoms are faithfully transferred 
to those of the emitted photons. The correspondence 
shown in Fig. |^ between the quasi-probability distribu- 
tion on the spin sphere and that on the complex-a plane 
indicates that a variety of photon states can be produced 
by merely rotating the spin vector of the SAS. We also 
found that this manipulation of few-photon states is pos- 
sible only if the atoms are in the SAS. Although the CAS 
can produce the photon-number squeezed state, the de- 
gree of squeezing is lower than that of the photon state 
produced by the SAS, and the phase squeezed state can 
never be produced by the CAS. The possible experimen- 
tal situations to implement our theory were discussed. By 
these schemes, we can generate the quantum-controlled 
few-photon state in the microcavity, and the quantum- 
controlled center-of-mass motion of trapped atoms. 
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APPENDIX A: VANISHING EXPECTATION 
VALUES IN THE JAYNES-CUMMINGS 
INTERACTION 



„iff* 



iEP' 



AH X ' 



\e~ *e>e 

i-iH' ot t 



It is assumed in Sec. IV that the expectation values 
(S x ) and (S2) always vanish if the time evolution is gov- 
erned by the Hamiltonian (Q) and when the initial state 
is I a) <£> \ S, M — S) with real a. In this appendix we give 
a general condition for this to be true. 

Since an expectation value of an Hermitian operator, 
say O, is real, it follows that 



iH Iot t\* 
/0 



0*e lH *) 



(Al) 



where the expectation values are taken with respect to 
the initial state |a) ® \S,A1 — S), and O* denotes an 
operator whose matrix elements are complex conjugates 
of those of O. In the second line of Eq. ( |Al[ ) we used the 
fact that the matrix element of the Hamiltonian (U2), 



(n\(S,M\H rot \S, M')\n') 



= gh y/K+ly/{S + M)(S-M + l)Sn, n '-iS M 



M' + l 



+V^V(S - M)(S + M+ l)<W+i<W'-i] . (A2) 

is real and hence H mt * = H mt . By a unitary trans- 
formation e z ^ Sz we have e ™s z jjrot e -nrS z _ _7jrot an( j 



e ivS '\a 



\a)®\S,M = S) 



a) ®\S,M = S). Applying 



this unitary transformation to the second line of (Al), 
the expectation value becomes 



t , a 



\ = (e iHVOtt e™ S *d*e 



iirSz 



-iH Ioz t\ 



0- 



(A3) 



Therefore, if e'^^'e"™ 5 ' = — <D, the expectation value 
( [A3] ) must vanish. The operators S x and &2 meet this 
condition. In general, an expectation value of an Hermi- 
tian operator that consists of operator products in which 
S x and a,2 appear odd- numbered times always vanishes. 
General conditions required for the initial state of the 
photon field pp and that of the atoms pa are p* F = pp 

and e mS "pAe 



PA- 



APPENDIX B: DERIVATION OF THE 



APPROXIMATE SOLUTIONS (35) 



In this appendix, we derive the approximate solutions 
([55|). It is convenient to define |3l| 



a. 



N 



1,2) 



~N 

t = gVNt, 



(Bla) 

(Bib) 
(Blc) 



in order to estimate errors of the approximation. The 
equations of motion for these operators have the forms 



d T S x 
d T S' y 

d T S' z 

drd^ 

d T a' 2 



-2a' 2 S' z , 

2{a' 1 S' y + a' 2 S' x ), 
-<?' 

s ~ . 



(B2a) 
(B2b) 

(B2c) 

(B2d) 

(B2e) 
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We assume that the initial state is the CAS \9 — 0, cj>) for 
the atoms and the coherent state |a) for the photon field, 
where a is taken to be real without loss of generality. 
Taking the expectation values of Eqs. (JB3) yields 



dr{S' y ) 

d T (S' z ) 
dr(a[) 



-2(a' 1 )(S' z )-2(Aa' 1 AS' z ) 
2(a' 1 )(S y )-2(Aa' 1 AS' y ), 

~(Sy), 



(B3a) 
(B3b) 
(B3c) 



where AO = O — (O) for any operators. It can be 
shown that (S' y ), (S' z ), and (a^) are of order unity, and 

(Aa'iA^) and (Aa^ASy) are of order l/N. If we neglect 
relative errors of l/N , the second terms of Eqs. (B3a) and 
(B3b) can be neglected, giving 



If wc set 



Or (S^ 

d T 0' z ) 

(S'y) 

(S' z ) 



2^ (^y), 



sin 9(t), 
cos 9(t), 



the equations of motion (|B3j) reduce to 
cP t 9{t) = sin 9(t), 



(B4a) 
(B4b) 
(B4c) 

(B5a) 
(B5b) 
(B5c) 

(B6) 



which has the same form as the equation of motion for the 
mechanical pendulum. The angular velocity of the pen- 
dulum corresponds to the field amplitude. The solutions 
for the initial condition 9(0) — and d T 9(0) — —2a', 
where a' = a/ \/~N , can be expressed in terms of Jacobi's 
elliptic functions as 



(S'y) 

(S' z ) 



—y/T—m sd(u\m)cd(u\m) , 
~ [2cd 2 Hm) - 1] , 

a nd(u\m), 
(S' x ) = 0, 



(B7a) 
(B7b) 

(B7c) 
(B7d) 



where u = t\1 + a' 2 and m = 1/(1 + a 12 ). 

The equations of motion for variances are written as 

d T ((Aa' 2 f) = -2(Aa' 2 AS' x ), (B8a) 
d T ((^S' x f) = -4(Aa' 2 AS' x )(S' z ) + 2(AS' x AS' z Aa' 2 ) 
+2(AS' z AS' x Aa' 2 ), (B8b) 
d T (Aa' 2 AS' x ) = -((A^) 2 ) - ((Aa' 2 ) 2 >(^) 

-((Aa' 2 ) 2 AS' z ). (B8c) 

It can be shown that the second-order fluctuations, such 
as ((AS' X ) 2 ), are of order l/N, and that the third-order 



fluctuations, such as (AS' x AS' z Aa' 2 ), are of order l/N 2 . 
Neglecting the third-order fluctuations in Eqs. (p38|), we 
have 



dr((Aa' 2 ) 2 ) 
dr((AS' x f) 
d T (Aa' 2 AS' x ) 



-2(Aa' 2 AS' x ), 
-4(Aa' 2 AS' x )(S' z 
-((AS' x f)~((A& 



(S' z ) 



(B9a) 
(B9b) 
(B9c) 



Using the form of (S' z ) in Eq. (B7b), which has at most 
a relative error of l/N, Eqs. (B8) reduce to the closed 
differential equations with relative errors 1/./V. They have 
three independent sets of solutions, and two of them are 
obtained as 



((AaM 2 



<(a4) 2 ) i = ! fsd^wa^u /»s 



i=nd 2 (zi|m) 



iV 



(Aa' 2 AS' x 



^5^sd(M|m)cd(u|m)nd(u|m) 



^nd 2 (u\m)E 2 (u\m) 



N 

^ [m sd(u\m)cd(u\m)E(u\m) + dn(u|m)] 

— [m sd(u\m)cd(u\m) E (u\m) +dn(it|m)]nd(ti|m) E (u\m) 

(BIO) 

The linear combination of these solutions to satisfy the 
initial conditions ((AS^) 2 ) = j^, ((Aa 2 ) 2 ) = and 



(Aa' 2 AS' x ) 



yields the solutions (35) 



APPENDIX C: A METHOD TO MINIMIZE THE 
PHASE FLUCTUATION 

In this appendix, we briefly show a method to ob- 
tain a photon state having the minimum phase fluctu- 
ation, which is the dotted curve in Fig. [l^. The variance 
of the Pegg-Barnett phase operator of the photon state 
J2 n Cn\n) F is given by 



{(Aor = — + 2 Y A ™,v 



(Cl) 



where A nm = (—l) n ~ m /(n — m) 2 . The coefficients that 
minimize the variance (Cl) satisfying the constraints 
c n = 1 an d J2n nc n = ™ are obtained by minimiz- 
ing the function 

F({c n }, A, 0) = 2 Y A nm c n c m + A I Y c 2 - 1 J 
+/3 (^Ynci-fiJ , 



(C2) 



where A and [3 are the Lagrange multipliers. The vari- 
ational problem dF/dc n = is equivalent to the eigen- 
value problem 



Y {2A nn > + nf36 nn i)c n/ + Ac„ = 0, 



(C3) 



1G 



which can be solved numerically. 
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